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Non-universal correlations due to localization are observed in statistical properties of experimental 
eigenfunctions of quantum chaotic and disordered microwave cavities. Varying energy E and mean 
free path I enable us to experimentally tune from localized to delocalized states. Large level-to-level 
Inverse Participation Ratio (IPR I2) fluctuations are observed for the disordered billiards, whose 
distribution is strongly asymmetric about {I2). The density auto-correlations of eigenfunctions are 
shown to decay exponentially and the decay lengths are experimentally determined. All the results 
are quantitatively consistent with calculations based upon nonlinear a— models. 
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The universal properties of quantum chaotic systems 
have been extensively studied in terms of their eigenvalue 
and eigenfunction statistics Q| . In Random Matrix The- 
ory (RMT), the Gaussian distribution of eigenfunction 
amplitudes ip(r) leads to the universal Porter-Thomas 
(PT) distribution for the densities |V>| 2 , which has been 
observed in microwave cavities pj as well as other sys- 
tems. However, non-universality has important manifes- 
tations, for instance due to periodic orbits which lead 
to scars in eigenfunctions, and localization arising from 
quantum interference in diffusion. While there have been 
many theoretical treatments, from semiclassical periodic 
orbit theories M to nonlinear sigma models ||, there 
have been few experimental studies of eigenfunctions be- 
cause of their lack of accessibility. 

In this paper, we present several striking manifesta- 
tions of localization in experimental eigenfunctions of dis- 
ordered microwave billiard cavities. Localization due to 
boundary or impurity scattering results in correlations 
that affect statistics and spatial correlations of the eigen- 
functions in several measures, leading to deviations from 
their universal values for quantum chaotic systems. The 
moments of the density distribution, /„ = J \fy(f)\ 2n d 3 r, 
in particular the Inverse Participation Ratio I2 (IPR), 
and their distributions Pj n (J n ), are important measures 
of the properties of the chaotic and disordered eigen- 
functions. In chaotic billiards, I2 has a mean value (I2) 
close to that of the universal 2-dimensional (2D) lim- 
iting value of 3.0, with small level-to- level fluctuations 
8I2 <C (^2)1 resulting in a nearly symmetric distribu- 
tion about (12)- In disordered billiards not only is the 
mean value (I2) ^> 3.0, but the fluctuations are also 
much greater SI2 ~ (h)- The IPR distribution for the 
disordered billiards is strongly asymmetric about (12), 
and is quantitatively consistent with the calculations 
based upon the nonlinear er— models of supersymmetry, 
parametrized by a conductivity g jU^. Spatial correla- 



tions are studied in terms of the density auto-correlation 
^| v E , (r)| 2 |^I'(r )| 2 ^ and are shown to die out more rapidly 
in the disordered billiards compared with the chaotic ge- 
ometries with a characteristic decay length given by the 
mean free path Here again the data are in quanti- 
tative agreement with the a— model calculations. Our 
results represent the first quantitative comparison of ex- 
periments and theory. 

The experiments were carried out using thin ( height d 
< 6 mm) cavities, whose cross-sections can be shaped 
in essentially arbitrary geometries. For these two di- 
mensional cavities, the operational wave equation is 
(V 2 + k 2 )^> = 0, where vb — E z is the microwave elec- 
tric field. Similar microwave experiments, which exploit 
this QM-EM mapping, have been used to study quan- 
tum chaos in closed HQ and open systems 0. Eigen- 
functions |'0( r )| 2 were directly measured using a cavity 
perturbation technique |J . Localization effects were ob- 
served by fabricating billiards in which 1cm circular tiles 
were placed in a 44cm x 21.8cm rectangle at random lo- 
cations (Fig.l) to act as hard scatterers. (The locations 
were generated using a random number generator and the 
tiles placed manually). Several realizations of the disor- 
dered geometry were experimentally studied, by varying 
the density of scatterers from 12 to 71, thus varying mean 
free path I ~ 8.9cm to 3.6cm. Earlier experiments 0] 
had shown that these disordered geometries are an ex- 
cellent experimental realization of the classic problem of 
an electron in a 2-D disordered potential. Subsequent to 
this work, there have been important theoretical develop- 
ments |^,|| , some motivated by our earlier experiments. 

The experimental eigenfunctions directly demonstrate 
the trend toward decreasing localization from disordered- 
chaotic-integrable as seen in Fig.0, which shows repre- 
sentative experimental eigenfunctions, along with their 
corresponding IPR I2, for several billiards. The strongly 
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localized state FigJ|(a) at / = 3.04GHz of the disordered 
billiard with N = 36 scatterers has very large I2 = 13.42. 
In contrast the more delocalized state at higher frequency 
/ = 7.33GHz (Fig.^(b)) explores almost all the available 
coordinate space similar to chaotic cavities and has a 
smaller I2 — 4.06. For the chaotic Sinai stadium, typi- 
cal values of the IPR are around 3.0 (I2 — 3.01 for this 
eigenfunction Fig.|l|(c)) while for the integrable rectangle 
billiard Figj^(d) IPR I 2 = 2.25 for all eigenfunctions. 
Fig.|l| demonstrates a key advantage of our experiments, 
which is that by varying I and wavevector k we are able 
to access a wide range of the disorder strength kl, from 
strongly localized states for kl < 1 to delocalized states 
with kl > 1. 

Eigenfunctions such as in Fig.|l| were then analyzed in 
terms of I2 and Pi 2 (12)- In the following, for convenience, 
we use the notation I2— J ^(r)! 4 ^ = J u 2 dv, w = (1% — 
3)/6 , u = |*(r)| 2 , and h= f\V(r)\ 2 dv = 1, and the 
integral is over the volume v (area in 2D ). Nearly 250 
eigenfunctions were analyzed each containing more than 
3200 eigenfunction values. 

Fig.|| shows level-to-level variations 12(E) for the Sinai 
stadium. Here I2 is mainly clustered around a mean value 
of (I2) = 3.0, with relatively small level-to-level fluctua- 
tions. The IPR distribution Pf 2 (/ 2 ) of the chaotic Sinai- 
stadium billiard is shown in Fig.||(top). Pi 2 {l2) is seen 
to be a nearly symmetric distribution with a small width 

h~{h) « (h)- 

Ref. Q| demonstrated that the Sinai-stadium billiard 
data obey the universal Porter-Thomas (P-T) distribu- 
tion P u (u) = (2nuy 1/2 cxp(-u/2), with u = j*| 2 to re- 
markable degree, while deviations from P-T were demon- 
strated due to localization in the disordered billiards. 
The IPR for P-T distribution can be immediately ob- 
tained I2 = J Q u 2 P u (u)du = 3.0, which is a universal 
value. Note that there are no fluctuations expected about 
this universal value in RMT, i.e., Pi 2 (l2) is a (5-function 
at I2 = 3 jio) . Clearly boundary scattering on the system 
length scale R leads to non-universal correlations (e.g. 
from periodic orbits leading to scars in the wavefunc- 
tions). This breaks the assumption in RMT of Gaussian 
fluctuations of the eigenfunction amplitude, and in turn 
leads to fluctuations in the distribution Pj 2 (I2) (although 
of narrow width) observed in Fig.^|(top) . 

Even more strikingly, the IPRs of the disordered bil- 
liards shown in Fig.^J display a remarkably large spectrum 
of level-to- level fluctuations (Fig.||). For small /, when 
A > I, strong localization leads to large values of I2 in 
the disordered cavity, which can be as high as ~ 20. It is 
worth noting that the density distributions P u (u) of the 
eigenfunctions deviate strongly from the P-T distribution 
and are consistent with the large IPR values. (In this pa- 
per we have focussed on the billiard with I = 5.1cm., but 
similar results were obtained for other billiards and will 
be presented in a larger publication.) 



As / is increased (or A is decreased), the IPR pro- 
gressively decreases, approaching the universal limiting 
value of 3, as shown in Fig.B. The corresponding distri- 
bution Pj 2 (l2) shown in Fig^bottom) is strikingly dif- 
ferent from the Sinai-stadium case, in that, it is strongly 
asymmetric, reflecting the very large I2 ^> (I 2 ) values 
that are present, and is strongly non-Gaussian. The IPR 
fluctuations in Fig.|^ are closely similar to the famous uni- 
versal conductance fluctuations of a mesoscopic system. 

Recently several theoretical studies have calculated the 
IPR distribution based on the supersymmetry method 
pLpJ. For a mesoscopic system the IPR distribution de- 
pends on the conductivity g of the system, defined as 
g = \n(R/l) I (w) , where R is the system size and I is the 
mean free path and (..) is the realization average for a 
fixed "disorder strength" 2kl. The resulting probability 
distribution P(I 2 ) for I 2 < (I 2 ) = 3 is f|: 

PiAh) = Cx|exp[-|(7 2 - <J 2 » - I c -f <*»-<*»»], (1) 

and the corresponding distribution for 7 2 3> (^2) '■ 

P l2 (h)=C 2 ^j-eM~9h) (2) 

where C\ and C 2 are normalization constants. 

The solid line in Fig.||(bottom) represents Eq.(2) for 
I2 ^> (I2) using a small conductivity g — 1.0 and is seen 
to describe the data very well. Another way of presenting 
the comparison with Eq.(2) is by using a scaled variable 
q = gl2 whence the distribution changes to a Porter- 
Thomas distribution in q. The system is within the weak 
disorder limit, i.e. the random phase approximation still 
valid, and g depends weakly on the disorder factor 2kl 
Jil"[ . Now rescaling the I 2 with g we obtain the distribu- 
tion P q (q) (unnormalized) and plot ln[P q (q)y/q] vs q in 
Fig.||. This shows a straight line which implies a good 
agreement of theory and experiment for IPR distribution 
in a moderate disordered region where J 2 3> (I 2 ). This 
is the first experimental observation of the asymmetric 
distribution predicted by Prigodin and Altshuler 0j. 

We now return to the case of the Sinai-stadium. Al- 
though a formulation in terms of a ballistic sigma model 
has been presented for chaotic cavities |[2| , it is not sim- 
ply amenable to experimental comparison. Instead we 
use Eq.l with the assumption that since I « R, a suitably 
large conductivity (g » 1) can be used. In our experi- 
mental case, g ~ 7.8 matches Eq.l not only for 7 2 < (^2), 
but also for I 2 > (^2), as shown in Fig. 3 (top). The nearly 
symmetric distribution can be understood since the fluc- 
tuations arise from correlations at the scale of the system 
length R, which is the only length scale in the problem. 
The values of g corresponding to the disordered billiards 
and the Sinai-stadium billiard are consistent with the (w) 
values and the mean free paths. 

We now show that the experimental data also di- 
rectly demonstrate the decay of spatial correlations 
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(^' 2 (r)\E' 2 (r')) of the eigenfunctions and obtain the de- 
cay length, which corresponds to the classical mean free 
path. To calculate the correlation with arbitrary disorder 
strength 2kl, defining K(r) — | Im G(r')\ 2 / (irv) 2 , where 
G(\r — r |) =< r\(E — i/) _1 |r' > is the Green function of 
the disordered system Hamiltonian, then it can be shown 
H that K(r) = |i j^Mkr(l + ^y)]dy\\ 

For the chaotic or ballistic limit (2kl >> 1), the result 
for Gaussian fluctuations is (ty 2 (r)^ 2 (?•')) ~ 1 + (1% — 
l)Jg(r — r'), where Jo is the first order Bessel function 
. The correlation for a moderate disordered case with 
correct limits can be derived repeating the calculations 
of Ref: ill 



(* 2 (r)* 2 (r')) ~ 1 + (J 2 - l)K(k\r - r'\) 
^l + (I 2 -l)J 2 (k\r-r'\)t 



(3) 
(4) 



These are valid when r — r' ^ /. We have solved 
K{k\r — r |) numerically. In Fig]5|, we plot the aver- 
age correlation derived from experimental data, numer- 
ical calculations of Eq.(3), and analytical expression of 
Eq.(4) for different disorder strengths 2kl. For the Sinai 
billiards, the average correlation starts at 3 and tend 
to 1 via Friedel oscillations, consistent with Eq.(4) with 
2kl = 37, which is very large and hence the result is close 
to that of the universal dependence given by kl = oo. 
For disordered billiards, the auto-correlation is very large 
(~ h) for short lengths, i.e. around \r — r'\ w due to 
localization, but the auto-correlation decays with a de- 
cay length scale, responsible for fast fall, and it should go 
to zero at large k\r — r'\. The experimental data shows 
good agreement with the numerical and analytical calcu- 
lations for moderate disorder, as shown in Fig.|B| for values 
2kl = 12 and 7. These values are in excellent agreement 
with the corresponding mean free paths (I — 3.6, 5.1, 
5.9 and 8.9) obtained by directly considering the num- 
ber of scatterers N, so that I ~ y/ ab/N and a and b are 
the sides of the enclosing rectangle. Thus our analysis di- 
rectly demonstrates localization and yields a quantitative 
measure in terms of the correlation decay length. 

Returning to the level-to-level I2 data in Fig.|| we 
note that they can be viewed as a form of localization- 
delocalization transition, tuned by the frequency / ! Rep- 
resenting the IPR as h(E) — h,sm{E) + SI 2 (E), i.e. a 
smooth part l2,sm(E) plus fluctuations 612(E), we dis- 
cuss the trends in 12(E). Calculations in Ref. based 
on infinite dimensional tight binding model shows that 
12(E) diverges exponentially near the critical point E c : 
h(E) = h(E = oc)exp(A/\E- E^), with \i = \ , due 
to very strong correlations of the wave function near E c . 
I2(E = 00) will be obviously the asymptotic universal 
value 3, as predicted by RMT in 2D. Our experimental 
IPR data is as large as IPR ~ 20 (strongly localized), 
decaying to ~ 4 (weakly localized) at high frequencies 
upto 10GHz. In the present case, the smallest scale of 



the system ~ 5cm sets a natural lower cut-off frequency 
fc = c/2l = 3GHz, so that there are no eigenstates for 
E < E c = f 2 . For E ^> E c expanding the above equation 
in first order, we obtain l2, S m(E) ~ 3 + B/\E — E c \^. 
In Fig.|| we have plotted this expression with /i = 0.5 
and B = 9.0 (obtained by a best fit) and compared 
with the experimental data. The above expression cap- 
tures the trend of the data. While an exact compar- 
ison with any expression is difficult since the fluctua- 
tions 512(E) > l2, S m(E), the analysis illustrates that we 
are observing a frequency driven localized to de-localized 
transition in a disordered medium in terms of the IPR. 

We have shown that the IPR is an extremely valu- 
able parameter to study real-space localization in quan- 
tum eigenfunctions. We have demonstrated for the first 
time a quantitative analysis of features well beyond uni- 
versality due to localization in experimental eigenfunc- 
tions. The observed IPR distribution is strongly non- 
Gaussian due to the correlations induced by scattering. 
The nonlinear sigma-model is in quantitative agreement 
with the experimental data for moderate localization. 
Our work thus provide experimental support for the ap- 
proach based upon quantum diffusion in the localization 
regime. 
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FIG. 1. (a,b) Experimental eigenfunctions of a disordered 
billiard with n = 36 scatterers (noted by the black dots) (a) 
a strongly localized state, / = 3.04GHz, I2 = 13.42 and (b) 
a delocalized state / = 7.33GHz, I2 = 4.06. Eigenfunctions 
of (c) Sinai-stadium billiard I2 = 3.01, (d) an integrable rect- 
angle (I2 — 2.25 for all states). 

FIG. 2. Large level-to-level fluctuations are observed in the 
IPR h{E) v.s. E of the disordered billiard. Also note the 
gradual trend towards a universal limiting value of 3.0 indi- 
cated by the solid line which represents a model described 
in the text. The fluctuations of the Sinai-stadium are much 
smaller and are clustered around < I2 >= 3.0. 

FIG. 3. IPR distribution P/ 2 (72) of the disordered billiard 
(bottom) is strongly asymmetric and non-Gaussian. The dis- 
tribution for the chaotic Sinai-stadium billiard (top) is nearly 
symmetric about the mean value 3.0. The lines represent cal- 
culations based on the nonlinear sigma-model. 

FIG. 4. The IPR distribution of Fig.3 replotted us- 
ing the scaled variable q = gh ■ The solid line is the 
Porter-Thomas distribution ln[P q (q)q 1 ^ 2 ] vs q and agrees for 
large h »< h >■ 



FIG. 5. Density auto-correlation (<l/ 2 (r)\E' 2 (r')) of eigen- 
functions of the Sinai-stadium and disordered billiards with 
fixed disordered strength 2kl. Experiment ( dotted lines), 
Eq.(3) (dashed lines), and Eq.(4) with 2kl =7, 12 and 37 
(solid lines). 
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